Induced Interactions for Ultracold Fermi Gases in Optical Lattices 
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We investigate the effect of optical lattices on the BCS superfluidity by using the Gorkov-Melik- 
Barkhudarov (GMB) correction for a two-component Fermi gas. We find that the suppression of 
the order parameter is strongly enhanced by the lattice effects. The predictions made by the GMB 
corrections are in qualitative and, for the cases studied, in quantitative agreement with previous 
quantum Monte Carlo results. We discuss how the GMB correction extends the validity of the 
mean- field theory to a wider range of tunable optical lattice systems in different dimensions. 
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Bardeen, Cooper, and Schrieffer (BCS) explained su- 
perconductivity by the condensation of fermion pairs 
in the presence of arbitrarily attractive interaction [![. 
Based on the idea of the BCS pairing, in dilute Fermi 
gases, the critical temperature was derived in terms of 
scattering length 0, IH . Gorkov and Melik-Barkhudarov 
(GMB) extended this calculation by incorporating many- 
body effects, which turned out to reduce the critical tem- 
perature by a factor (4c) 1 / 4 « 2.22 0]. Fermionic super- 
fluidity has recently attracted renewed attention in con- 
nection with the realization of ultra-cold atomic gases 
that allow direct observation of quantum many-body 
phenomena in highly controllable environments 0, @, 0] • 
In particular, optical lattices are a perfect platform 
for emulating crystalline structures of superconductors. 
While indirect evidence of superfluidity in a system with 
an optical lattices has been recently reported [8| j the full 
characterization of fermionic superfluidity and strongly 
correlated quantum states in optical lattices is still under 
active study, both theoretical and experimental [Io| . 
In this Letter we focus on how the lattice potential influ- 
ences the BCS-type superfluid transition by employing 
the GMB correction. 

We calculate the mean-field BCS order parameters at 
zero temperature in three and two dimensional (3D and 
2D) lattices with various settings including the crossover 
from 3D to ID. In all the ranges of lattice parameters ex- 
amined, we find that the induced interaction introduced 
by the correction leads to remarkable reduction in the 
order parameter from the usual BCS result. This de- 
viation turns out to be much more pronounced in the 
lattices than in homogeneous gases and becomes increas- 
ingly significant at higher fillings. In particular, in 2D, 
we find quantitative agreement with previous quantum 
Monte Carlo (QMC) calculations for the cases studied. 
Furthermore, near half filling in 2D, the rapid decreasing 
behavior of the order parameter is in qualitative agree- 
ment with the QMC predictions. At half filling in 2D, 
the induced interaction diverges because of Fermi surface 
nesting. This divergence is connected to the signature of 
the charge density waves, known to coexist with super- 



fluidity at half filling in 2D. 

We consider a system composed of two different 
fermionic species denoted by f and |. Each compo- 
nent is in a lattice with adjustable tunneling strengths 
ty a and t| Q in direction a G {x, y, z}. When the lat- 
tice potential is sufficiently deep so that we can consider 
only nearest-neighbor tunnelings and on-site interactions, 
the system is described by the Hubbard Hamiltonian 

H = - E ff , a J2 la tvaK-aic + U J^i "Ti^ii ~ A 1 Ea,i ^ > 

where K ai = i>\ i+1 i>ai + ^Vvi+i, and is the 

annihilation(creation) operator for atoms of type a at a 
site i = (ix,iy,iz)- The chemical potential and density 
operator are denoted by [i and n CT) i, respectively. We 
consider negative interaction strengths f7o- For nonin- 
teracting gases, the above Hamiltonian is diagonalized 
with dispersion £ CT (k) = 2 ^2 a t aa [1 — cos(fc a )]— fi, where 
the lattice spacing is chosen to be unity. In the weak 
coupling regime, using the standard BCS theory with 
this dispersion and the interaction U$, without many- 
body corrections, one recovers the usual BCS predic- 
tion for the critical temperature in the long wavelength 
limit. However, because the scattering length a gives 
exponential contribution to the critical temperature as 
T c oc exp(— TT/kp\a\), where kp is the Fermi momen- 
tum, a small correction to the interaction term fcpa 
can considerably change T c even in the weak coupling 
regime. For instance, a second-order correction 5 in 
k F \a\ -> k F \a\ (1 + 8k F \a\) leads to T c -> e s T c . 

For a two component Fermi gas with an s-wave in- 
teraction between components, the relevant second order 
correction to the effective interaction is represented by 
the diagram in Fig. [IJa), which describes the exchange 



of density and spin fluctuations [111 . |12[ . The diagram 
leads to the induced interaction term, which can be de- 
rived for an infinite-size system in D dimensions as 



C/ ind (p,k) = -ul 



rfq 



/t,p+k+q - /j.,, 



(2^£T(P + k + q)-£i(q)' 



(1) 



where the Fermi distribution f a ^ = 1/ [1 + exp(/3£ CT (k))], 
with p = l/fc^T. There are two noticeable properties in 
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FIG. 1: (a) The diagram representing the induced interaction 
f/in<j(p, k). Arrowed and dashed lines describe fermionic prop- 
agators and the coupling Uo between the atoms, (b) Equiv- 
alence between the zero-temperature order parameter A and 
the critical temperature T c in the weak coupling regime. 



Eq. |T]). First, U[ n d is always positive. Thus, this cor- 
rection screens the negative interatomic potential, which 
consequently reduces the critical temperature. Second, 
the static Lindhard function, representing the spin or 
density susceptibility Xo of a noninteracting gas, is found 
in U- m d if there is no difference in the energy dispersion 
between the components, or simply t^ a — t^ a , which is 
the case that we mainly consider here. 

In the weak coupling regime, the induced interac- 
tion correction near the Fermi surface dominantly con- 
tributes to the calculations of the BCS order parame- 
ter and the critical temperature (ill . [l2j ]. and then the 
effective interaction is approximately given by only the 
Fermi surface momenta. Averaging f7i n d(p,k) over the 
Fermi surface, the induced interaction becomes (U-md) = 
|g j^g | fg dS p J s ^ dSkUi n d(p, k), where S a denotes the 
Fermi surface of the component a and 1 5*0- 1 = J s dS is 
the area. Finally, the effective interaction is written as 
U c tf = Uo + (Uind)- This effective interaction replaces 
the interatomic interaction in calculations. Having other 
parameters fixed, U e g becomes zero when Uo equals to 
U c = — C/o/(E/i n d}- The susceptibility of an interacting 
gas is given as Xo/(l + UqXo) within the random phase 
approximation, where xo is the susceptibility of the non- 
interacting gas. At Uo = U c , this quantity diverges, of- 
ten indicating a possibility of charge ordered phase. Our 
approach is formally valid below U c , and with the GMB 
correction, the criterion for the weak coupling regime can 
be established by U c <C Uq < 0. 

We calculate zero-temperature order parameters by us- 
ing this effective interaction U c g in the mean-field for- 
14j . The diagram in Fig. QJa) representing 
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the induced interaction is relevant in normal states, and 
thus the critical temperature T c has been of interest in 
previous studies. However, in the weak coupling limit, 
the validity of the induced interaction can be readily ex- 
tended for the calculation of the zero-temperature order 
parameter A because the contribution of the broken sym- 
metry phase is expected to be of higher order in A that 
becomes very small. Moreover, with U = Ueff, we con- 
firm that a general relation 2A/kBT c ~ 3.53 [H,[I3| still 
holds in lattices, as shown in Fig. [Hb) , in the weak cou- 
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FIG. 2: Three-dimensional lattices. The order parameters 
A (with the GMB correction) and A (0) (the usual BCS re- 
sult) are calculated with Uo = —3t as a function of chemical 
potential p>. The comparison between A and A^ ' shows an in- 
creasing deviation as [i increases. The ratio A^°'/A becomes 
~ 25 at half filling (fi = 6t). The inset shows comparison 
between the magnitude of the interactions. Here and in the 
other figures, the critical coupling \U C \ (see text) allows to 
estimate the ranges of \Uo\ and fi where the BCS mean-field 
theory with the GMB correction is applicable. 



pling regime. The calculation of A is computationally 
less demanding than T c but gives a concise picture of 
the transition. At very low filling factors (small /i), one 
recovers the well-known prefactor 2.22 in comparison be- 
tween the usual BCS result and our calculation with the 
correction in f7 e ff- I n contrast, as the filling factor in- 
creases, the reduction of the order parameter with the 
correction becomes much more subtle in a lattice than 
indicated by the prefactor in continuum. 

Figure [5] shows the effect of the induced interaction 
in isotropic 3D lattices. We find that the order param- 
eter A with the correction shows a dramatic deviation 
from the usual BCS result A^ without the correction. 
The ratio A< ) /A turns out to be nearly 25 at half filling, 
which implies that the usual BCS prediction largely over- 
estimates the order parameter and critical temperature. 
While the order parameter A with the correction is max- 
imized around /i = At where the Fermi surface reaches 
the Brillouin zone boundaries, it decreases substantially 
at higher filling factors because of increasing contribu- 
tions of states with p + k + q in Eq. |T|) outside the 
first Brillouin zone. This deviation from the usual BCS 
result that we find here is qualitatively consistent with 
the previous study of the 1/D correction in high dimen- 
sions 



171 ] , where the order-of-unity reduction of the or- 
der parameter was estimated at half filling. However, our 
calculations in 3D lattices reveal much more significant 
suppression in A with the GMB correction. In all the 
figures, the curve for |J7 C | provides a simple estimate of 
the area in the parameter space (Uo, /i) where the GMB 
correction can extend the validity of the BCS mean-field 
theory. The range of Uo in this space decreases as (i in- 
creases towards half filling in 3D and 2D lattices, which 
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may affect the accuracy. Particularly in 2D, half filling 
is not included in this space because of divergent (t/; n d)- 

In 2D lattices, the order parameter A decreases very 
rapidly near half filling as the screening by the induced 
interaction dominates (see Fig. [3]). This rapid decrease 
of A near half filling is in agreement with the pre vious 
QMC results on the critical temperatures [l8L Il9|. In 
contrast, the usual BCS mean-field calculation without 
the correction suggests a monotonically increasing order 
parameter A^ ) when approaching half filling. At lower 
filling factors, the order parameter with the GMB cor- 
rection turns out to be around 5 times smaller than the 
one without the correction, which is consistent with other 
previous estimations of many-body effects [IH . At 
higher fi close to half filling, the deviation from the usual 
BCS theory becomes even more substantial and leads to 
highly suppressed order parameter near half filling. 

At half filling, it is known that the 2D attractive Hub- 
bard model has the charge-density-wave order and the 
pairing order coexisting in the ground state, and the crit- 
ical temperature of the superfluid transition goes to zero. 
Our calculations show that the induced interaction loga- 
rithmically diverges because of the Fermi surface nesting, 
making the denominator in Eq. ((T|) infinitesimally small 
with the nesting vector p + k = (±7r, ±7r), for all q at the 
Fermi surface, by mapping one side of the Fermi surface 
onto the other side [see Fig. [3f d)] . Despite the tendency 
for the order parameter A to vanish when approaching 
half filling, the order parameter is not well defined with 
the large correction in our perturbative approach. How- 
ever, this divergence of the induced interaction means 
that, for arbitrary small Uq, the susceptibility diverges. It 
thus provides the connection to a different type of phase 
that cannot be anticipated by the usual BCS mean-field 
theory. The divergent Lindhard function appearing in 
the correction term can be interpreted as the signature 
of the charge density waves [12] ■ 

For direct comparison with the QMC values of T c in 2D 
lattices, we have used Uq = —44. With the GMB correc- 
tions, we obtained A (T c ) ~ 0.074 (0.044), 0.034 (0.024), 
0.008* (0.005t) at the filling factors (|t/ c |'s) of 0.2 (7.34), 
0.25 (6.24), 0.3 (5.44). The farther away U c is from 
Uq, the more accurate A is expected. The QMC result 
IsbTc ~ 0.054 at quarter filling 18, [l9| is remarkably 
close to our value of A, and the results in are of the 
same order of magnitude as ours. Note that the usual 
BCS mean-field calculations for these parameters would 
give results that are about 10-35 times larger than the 
QMC and our GMB corrected mean-field results. 

Motivated by the fact that anisotropy is easily con- 
trollable in optical lattices, we now explore dimensional 
crossover from 3D to ID by introducing directional dif- 
ference in the tunneling strengths t a . For this purpose, 
we define the lattice anisotropy as a ratio of the tunneling 
strengths, 4 = t y /t x = t z /t x , with which one can change 
the dimensionality from 3D (4 = 1) to ID (4 = 0). Fig- 
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FIG. 3: Two-dimensional lattices, (a) The effective interac- 
tion U e ff and (b) the order parameter A are calculated with 
Uo = —1.54. The usual BCS mean-field result A' ' shows very 
large deviation from A with the GMB correction, exhibiting 
A (0) /A « 10 near half filling n = 44. (c) Logarithmic di- 
vergence of the GMB correction and (d) Fermi seas (shaded 
area) at half filling indicating the nesting of Fermi surfaces 
(dashed lines) with momentum transfers p + k = (±7r, ±7r). 



ure [4] shows the effect of the lattice anisotropy on the 
induced interaction and the order parameter. As the 
anisotropy evolves with 4, the screening by the induced 
interaction becomes stronger and finally diverges in the 
limit of ID because of Fermi surface nesting. 

Similar to isotropic 3D cases, the order parameter A 
in the anisotropic lattice also shows highly suppressed 
values compared with the usual BCS result A^ ^. While 
the transition from 3D to ID appears continuous in the 
induced interaction (C4 ind ), we identify two special points 
of 4 indicating structural changes of the Fermi surface. 
First, (t/ind) has a kink around 4 = 0.5 at which A begins 
to decrease. The Fermi surface is closed originally in the 
3D lattice with the given chemical potential \i = 2t x . 
With decreasing 4, the Fermi surface becomes deformed, 
and then at 4 = 0.5, the Fermi surface becomes open. 
The second is a bump of (t/ind) near 4 = 0.25 at which a 
dimensional change of the Fermi surface occurs and the 
nesting effect develops to escalate (tAnd)- As plotted in 
Fig.HJc), finally at 4 = 0.25, the surface completely opens 
in k y and k z directions and splits into two disconnected 
sheets causing the nesting effect. 

However, in the quasi-lD regime, it turns out that the 
parameter space given by U c does not cover the low 4 re- 
gion, and our calculation predicts a vanishing order pa- 
rameter at low 4, which deviates from previous rigorous 
studies of the Hubbard model. In the attractive Hubbard 
model in quasi-lD, the spin gap and the critical tempera- 
ture are finite [3]. In the limit of ID, the gap is still finite 
though the critical temperature goes to zero 25j . Singlet 
superfluidity dominates in the ground state, but there is 
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FIG. 4: Crossover from 3D to ID. (a) The induced interac- 
tion (t/ind) and (b) the order parameter A as a function of 
the lattice anisotropy t = t y /t x — t z /t x in anisotropic three- 
dimensional lattices. In (b), the usual BCS mean-field results 
A' ' without the correction is given for comparison, (c) Fermi 
surfaces projected to k x -k y space at t = 0.26, 0.25, 0.24, 0.05 
(from center). The chemical potential is fixed at /i — 2t x and 
Uo = — 3t x is used in the calculations. 



ings suggest that the effective theory with the many-body 
correction to the interatomic interaction can significantly 
extend the applicability of the mean-field calculations in 
the lower dimensions, namely 3D and 2D, and in the 
crossover from 3D to ID lattices, in spite of the obvious 
failure in the strict ID limit. With the GMB correction, 
the simple mean-field calculation can also provide quan- 
titatively reliable values in a wider range of the coupling 
strength, without sophisticated QMC calculations. The 
unanticipated large suppression of the order parameter 
at high filling factors highlights the practical importance 
of our results, which may provide a new insight to the 
issue of the critical temperature in future realizations of 
fermionic superfluids in optical lattices. 

The authors thank Dr. T.K. Koponen for fruitful dis- 
cussions. This work was supported by Academy of Fin- 
land and EuroQUAM/FerMix (Projects No. 213362, 
217041, 217043, and 210953) and conducted as a part 
of a EURYI scheme grant jig. 



no true long range order in the ID Hubbard model [261 ]. 

We have also considered the problem of the fermions in 
component-dependent lattice potentials [27J where each 
component experiences a different tunneling strength in 
a lattice. This difference in tunneling in a lattice is anal- 
ogous to unequal effective masses of Fermi gases in con- 
tinuum. In 3D lattices, we have found that the screen- 
ing effect of the induced interaction becomes stronger 
as the difference between the tunneling strengths in- 
creases, which agrees well with the results for homoge- 
neous gases 
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28] where similarly the stronger screen- 
ing effect at the larger mass imbalance was found. 

In conclusion, we have found that the presence of the 
optical lattices substantially strengthens the effect of the 
GMB correction on the BCS superfluidity. The conse- 
quent suppression of the order parameter is found to be 
much beyond the ratio 2.22 predicted in homogeneous 
gases, which agrees with the estimations of other previ- 
ous many-body correction studies at low filling factors. 
As the filling factor becomes higher, the inclusion of the 
correction becomes increasingly important. For instance, 
the order parameter turns out to be almost 25 times 
smaller with the correction than the usual BCS mean- 
field results at half filing in 3D lattices. Moreover, the 
behavior of the order parameter in 2D lattices shows ex- 
cellent agreement with the previous QMC values. Natu- 
rally, when the correction becomes very large, our pertur- 
bative approach breaks. The divergence of the correction 
is related to the phase at half filling in 2D where super- 
fluid order and charge-density- wave order coexist. 

One of the general shortcomings of a mean-field the- 
ory is that it gives a valid approximation only in high 
dimensions. Particularly for BCS superfluidity, our find- 
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